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Abstract 

This work is devoted to the development of a novel theoretical approach, 
named hybrid approach, to handle a localized bottleneck in a symmetrically 
conpled two-channel totally asymmetric simple exclusion process with Lang- 
mnir kinetics. The hybrid approach is combined with singnlar pertnrbation 
techniqne to get steady-state phase diagrams and density profiles. We have 
thoronghly examined the role played by the strength of bottleneck, binding 
constant and lane-changing rate in the system dynamics. The appearances 
of bottleneck-induced shock, a bottleneck phase and Meissner phase are ex¬ 
plained. Further, the critical values of bottleneck rate are identified, which 
signify the changes in the topology of phase diagram. It is also fonnd that 
increase in lane-changing rate as well as nneqnal attachment, detachment 
rates weaken the bottleneck effect. Our theoretical arguments are in good 
agreement with extensively performed Monte Carlo simnlations. 
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1. Introduction 


In protein synthesis, the genetic information is deciphered into proteins 
by molecular machines called ribosomes that attach themselves at the start 
end of mRNA, move aloi^ the chain in a unidirectional manner and finally 


detach at the stop end 


2j. Each translation step requires the binding of 


a freely diffusing trans; 
specihc to each codon 


er-RNA (tRNA) molecule, carrying the amino acid 


,Q. 


The important factor affecting the ribosome 
translation rate is relative concentrations of tRNA, which may vary from 
codon to codon. The codons with lower concentrations of tRNA reduce 


the protein synthesis ra 
a homogeneous system 


;e and 
32, 3, 


hus plays the role of an inhomogeneity in 


occur natural 
blood flow 


34j |. Apart from this, inhomogeneities also 
y in many other transport systems such as vehicular traffic 11| . 


24j | and fiow of data in a Von Neumann architecture 20(]. In traffic 


ffow, the ongoing construction on roads, a slow moving vehicle or an accident 
can lead to slow down the flow rate on highways and can lead to congestion. 
Further, the separation of the CPU and the memory in computers creates 
Von Neumann bottleneck, which limits the performance of the computer via 
limited bandwidth between the CPU and the memory. 


Totally asymmetric simple exclusion process (TASEP) 3|, l2l| is well 


known to be a paradigmatic model for studying stochastic transport in 
many-particle systems. Both single-channel and multi-channel TASEPs have 


been well e xp 

a lia y, bs 


ored theoretically as well as using Monte Carlo simulations 


fl 


31| . Further, one has to take into account the fact that 


the proteins as molecular motors can also attach from the bulk reservoir 


or detach from it (Langmuir Kinetics (LK)) 


Q- 


In contrast to TASEP, 
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where total number of particles remain conserved, the additional attachment- 
detachment dynamics (LK dynamics) violate the particle conservation and 


leads to many interesting phenomena 


In literature, the consequences 


of coupling of the two different dynamics : TASEP and LK have been well 


analyzed in single-channel 


22 


23 


systems 


25| as well as two-channel homogeneous 


I lid, lid, l3d] • The idea of slowing down of particles at certain defect 
positions can be incorporated in the form of a set of inhomogeneous lattice 
sites (bottleneck), either as a single unit or randomly distributed over the 
whole lattice in TASEP. This type of disorder is known as site-wise disor¬ 
der. Another type of disorder studied in literature is particle-wise, where 
a slow moving particle itself acts as an inhomogeneity in the system. The 
present work focuses on site-wise disorder which is suitable to model the 
inhomogeneities present in transport of molecular motors. 

Although a lot of work has been done on homogeneous TASEPs, the effect 
of disorder on the steady-state dynamics of such systems is not well under¬ 


stood. Several studies have been performed on sin 


TASEPs with 


26 


29| as well as without LK 


e-channel in 
7, Is, 14, 15 


homogeneous 


18 


While investigating the role of a bottleneck in a closed TASEP 


,l35|. 


it was found that even the presence of a single bottleneck site can pro¬ 
duce shock profile and a plat eau in the fundamental current-density rela¬ 


tion. Later, Kolomeisky 


18l | examined even richer case of single-channel 


open system and explored the consequences of the inhomogeneity with both 
faster and slower transition rates. He divided the system into two homoge¬ 
neous TASEPs coupled at the single bottleneck site (defect mean-field theory 
(DMFT)) and proved analytically that a fast site has no effect on the phase 
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diagram; whereas a slow site leads to shifting of the phase boundaries only. 
He also tested the theoretical results with Monte Carlo simulations and found 
a good agreement in low density (LD) and high density (HD) phases; while a 
little deviation in maximal current (MC) phase. Another analytical approach 
namely hnite-segment mean-held theory (FSMFT) was introduced by Chou 


and Laktaos 


3 to study clusters of slow codons in protein synthesis. They 


found that ribosome density prohles near neighboring clusters of slow codons 
suppress the proteins synthesis. Dong et ah ^ generalized the DMFT to 
study the effect of two bottlenecks on the protein production rate and also 
performed extensive Monte Carlo simulations to conclude that the location 
as well as spacing between the bottlenecks affect the production rate, point¬ 
ing out an important phenomenon namely edge effect. The investigation 
into edge effects was further carried out by Greulich and Scadshneider ^ 
to generate the phase diagrams of inhomogeneous TASEP using interacting 
subsystem approximation (ISA). They could successfully explain the interac¬ 
tions of defects with the boundaries of the single-channel system. The more 
complex case of inhomogeneous TASEP in the presence of Langmuir kinetics 
was studied by Qiu et ah j29|. They calculated phase diagrams and density 


prohles by adopting the concept of DMFT and also studied the effect of s 
hopping rate and detachment rate on the phase diagram. Pierobon et ah 


ow 


26| 


provided a detailed study on the role of a bottleneck in a TASEP with LK, 
using an effective mean-held theory and Monte Carlo simulations. They in¬ 
troduced the concept of carrying capacity to identify various novel phases 
called bottleneck phases. Importantly, all of the above studies focused on 
single-channel inhomogeneous systems. 
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Parallel to the inhomo^neous single-c 


channel transport systems |g, l9|, [lO, Il6l. 


28|, 


lannel systems, particles in mnlti- 


3l| may also confront a bottleneck, 


present in either one or in more than one channels. The importance of stndy- 
ing the mnlti-channel system lies in the fact that it can act as a framework 
for extending the analysis to networks. Dne to the complexity in dynamics 
generating from the interactions between different channels, it is difficnlt to 
examine the effects of inhomogeneity in a mnlti-channel open system. Up 
to onr knowl edg e, the only contribntion in this direction has been made by 


Wang et al. 


37l |. which explored the effect of a local inhomogeneity in one 


of the lanes of a two-lane TASEP with LK nnder a symmetric lane changing 
rnle. They extended the DMFT to a two-channel system by incorporating 
the concept of effective injection and removal rates at the inhomogeneons 
lattice site. Despite a good agreement between the solntion of the mean-held 
eqnations and Monte Carlo simnlations in ref. , this approach fails to pro- 
dnce the steady-state phase diagrams of the two-channel system. Moreover, 
it is not feasible to analyze the role played by varions parameters snch as 
lane-changing rate, attachment-detachment rate and strength of bottleneck 
on the steady-state dynamics. One can infer from here that thongh DMFT 
is capable to provide analytic solntions for single-channel inhomogeneons 


TASEPs with 


26 


29| as well as withont LK j^, it lacks some important in¬ 


gredients to generate a complete pictnre of the dynamics of the corresponding 
two-channel system as discussed in section [31 This motivates us to develop 
a new approach to handle the bottleneck in a two-channel TASEP with LK, 
which not only overcomes the existing limitations, but also demonstrates the 
unexplored dynamics of two-channel inhomogeneous systems. 
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The objective of the proposed study is two-fold, to develop a general 
theoretical approach, which is capable to produce the phase diagrams and 
to study the effect of various system parameters on the steady-state phases. 
In this paper, we attempt to provide a complete picture of the dynamics 
of two-channel symmetrically coupled TASEP with LK in the presence of 
a single localized bottleneck in one of the two channels by adopting a new 
and simplihed approach, called the hybrid approach. We have also validated 
the theoretical results with Monte Carlo simulations. The paper is organized 
as follows. In section |2l we dehne the model under examination and the 
governing dynamical rules. We briefly discuss the limitations of the earlier 
approaches in sec. [3l The theoretical hybrid approach and Monte Carlo 
simulations are covered in section H] and section |3 A thorough analysis 
of the stationary properties of the model is discussed in section [6l In the 
concluding section [TJ we summarize the results and future prospectives of 
our work. 

2. Two-channel inhomogeneous TASEP with LK 

We dehne our model in a two-channel (L, 2) lattice, where L is the length 
of a channel. The two channels are denoted by A and B, in which particles 
are distributed under hard-core exclusion principle (See Fig. [T]). We adopt 
random-sequential update rules for the dynamical evolution of the system. 

For each time step, a lattice site (i, j); i = 1, 2, 3,. L; j = A, B is randomly 

chosen. The state of the system is characterized by a set of occupation 

numbers Tij {i = 1, 2, 3,.L; j = A,B), each of which is either zero (vacant 

site) or one (occupied site). At entrance {i = 1), a particle can enter the 
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lattice with a rate a provided rij = 0; and at exit (i = L), a particle can 
leave the lattice with a rate (3 when tlj = 1. In the bnlk, if Tij = 1, then the 
particle at the site {i,j) hrstly tries to detach itself from the system with a 
rate (detachment rate) and if it fails then it moves forward to site (z + 1, j) 
with a rate Pij provided Tj+ij = 0; otherwise it attempts to shift to other 
lane with a rate uj, only if the target site is vacant. On the other hand, if 
Tjj = 0 ; i = 2, 3,L — 1, a particle can attach to the site {i,j) with a rate 
uja (attachment rate). Here, horizontal transition rate Pij is inhomogeneous 
and is given by the following binary distribution 



q ] i = m Sz j = A 
1 ; otherwise 


( 1 ) 


It is clear from Eq. CDi that the bottleneck connects the sites sX i = m 
and i = m + 1 in lane A (Fig. [1]). Here, q denotes the transition rate of a 
particle on passing through the bottleneck and will be called as bottleneck 
rate throughout this paper. We wish to analyze the effects of a localized 
bottleneck in the bulk with no boundary interactions of the open system. 
This adds up to the assumption viz., 1 << m << L. 

3. Inadequacies of DMFT 

The only step to handle a single-site bottleneck in a two-channel TASEP 

n 

with LK was put forth by Wang et ah [37[. They generalized the concept 
of DMFT from a single-channel inhomogeneous TASEP with LK to the cor¬ 
responding two-channel system. Although inhomogeneity is present in only 
one lane, yet they divided both the lanes into two sub-lanes separated at the 
position of bottleneck site and introduced effective exit and entrance rates 
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Lane B 


Lane A 


Figure 1: Schematic diagram of a symmetrically coupled two-channel TASEP with LK 
with bottleneck in lane A at site. Crossed arrows show forbidden transitions. 

at bottleneck site, respectively in each lane. Using mean-field approximation 
with some simplified assumptions, they computed the nonlinear system of 
equations in terms of effective rates and solved it numerically. The steady- 
state density profiles obtained using DMFT are found to be in good agree¬ 
ment with Monte Carlo simulations. But the incompetency of this approach 
lies in the fact that the nonlinear system may not have a unique feasible 
solution always for all choices of parameters, due to which the complete in¬ 
formation about the steady-state dynamics is missing. Moreover, the phase 
diagrams cannot be obtained and it is not feasible to quantify the role played 
by different parameters such as lane changing rate, attachment-detachment 
rates and bottleneck rate etc. Since this approach uses the concept of divid¬ 
ing each lane (inhomogeneous as well as homogeneous) into subsystems, its 
generalization either to multi-channel system or to include multiple defect 
sites would be a cumbersome job to handle. 
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4. Theoretical approach 


As discussed in the previous section that DMFT is incapable of explaining 
certain important aspects, we introduce a novel hybrid approach to deal with 
the bottleneck in a multi-channel TASEP. Instead of dividing both the lanes 
into two sublanes each, we divide only the inhomogeneous lane (lane A) into 
two subsystems viz., 1 <i < m and m + 1 < i < L, connected through sites 
i = m and i = m + 1. 

Firstly, we compute the temporal evolution of occupation probabilities 
in the bulk (1 < i < L, j = A, B) from the following set of master equations: 

- Aj)) - - A+ij)) + - Aj) - ^d{ri,j) 

T^(^Ti^ATi+l,AO- - B,b) - (A,BA+i,i3(l - 

( 2 ) 

where (• ■ ■) denotes the statistical average and last term on right-hand side 
takes a negative and a positive sign for lane A and B, respectively. At the 
boundaries, the particle densities evolve according to 

= a(l - Tij) - Pl,j{Tij{l - T2,j)), (3) 

= PL,j{rL-i,j{l - tlj)) - ( 4 ) 

Factorizing the correlations using mean-held approximation, we get (rjjTj+ij) 


4-1. Continuum mean-field equations 

To hnd the continuum limit of the master equations, we dehne lattice 
constant e = 1/L, rescale the time and other kinetic rates as t' = t/L, 


CLd = oOdL and = uL 


Id, 


25|. 
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Writing x = (i — l)e , we replace binary discrete variables Tij with contin¬ 
uous variables pij G [0,1] and retain the terms up to second-order in Taylor’s 
series expansion to obtain 


Pi±ij - Pij ± + 


^P^,j , ^ ^^^ 3 ^ 


(5) 


dx ' 2 dx'^ 

4-1.1. Continuum mean-field equation for lane B 

Since lane B is free from any inhomogeneity, we can drop the subscript i 
and compute the following time evolution equation from master equation ([2]) 
for average density in lane B, denoted by 


~ ~ Pa))j 


( 6 ) 


where K = Qa/^d is the binding constant and an important parameter in 
analysing the system dynamics 251. 


4 . 1 . 2 . Hybrid mean-field system of equations for lane A 

We divide lane A into two sublattices viz., 1 < i < m and m-|- 1 < i < L, 
connected through sites i = m and i = m + 1. The two homogeneous 
subsystems connected by the bottleneck in lane A will be referred to as 
left and right subsystem. Since, each of the two subsystems is individually 
a homogeneous TASEP with LK, we can End the continuum limit of the 
mean-held approximate equation for each subsystem in lane A by proceeding 
in a similar fashion as for lane B. Keeping the master equations at (m. A) 
and {m -\- 1, A) sites intact after rescaling of aforementioned variables, we 
obtain the following hybrid system of equations for time evolution of average 
density in lane A, denoted by pa- 


10 







Hybrid system 

Continuum part: For x G (0, me) lj((m + l)e, 1), we have 


dpA 

dt' 


^ + jriPA ~ pa) + ^d{K — {K + 1 )pa) — Vl[p\{1 — pb) 


( 7 ) 


( 8 ) 


2 dx"^ ' dx 
-PbC- pa))- 
Discrete part: 

"OP — ~Pm — l,AC Pm, a) ~Pm,AC Pm+1,A) ~t“ (i^ + l)Pm,A) 

^{,Pm,APm-\-l,AC Pm^s) Pm,BPm-\-l,BC Pm,A}) i 

“ ~Pm,AC Ptu+I-Al) “Pm+1,A(1 P?71+2 .a) “t“ (^K -j- \)PjyiJ^\ a) 

or e e 

^ (Pm+l,APm+2,A(1 Pm+im) Pm+l,_BPm+2,S (1 Pm-\-l,A)) • 

In the continuum limit, the boundary equations ([3]) and (jl]) reduce to Py4(0) = 
Ps(0) = a and Pa(1) = PsiX) = 1-/9. Along with these boundary condi¬ 
tions, the system of equations ([S]), dZ]) and ([HD is consistent for which, a 
unique feasible solution always exists. With this favour, our new approach 
overcomes the limitation of DMFT, which was unable to guarantee a unique 
feasible solution to the nonlinear system of equations 37|. Moreover, no ap¬ 


proximation, other than mean-field, has been used to generate the hv 


system, which is again an advantage over the methodology of DMFT 


37 


Drid 

for 


a two-channel inhomogeneous TASEP with LK. It can be easily seen that the 
hybrid approach only segments the inhomogeneous lane into two parts, with¬ 
out disturbing the homogeneous lane. This idea makes it applicable to more 
general systems such as multi-channel or networks and those with multiple 
bottleneck sites. 


4-2. Steady-state solution 

We introduce a new space variable x to merge the equations ( 0,0 and 
(ISP into a single system, which further reduces to the following system in 
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steady-state, 
e (PpA 


,dpA 


( 9 ) 


2 dP + - 1)^ + ^diK -iK + 1 )pa) - ^p\{l - Pb) + ^P%{1 - pa) = 0, 

—1,a(1 Pm, a) QPm,A{^ Pm+1,A) “t “t ^)Pm,A) 

^Pm,APm+l,.4{l Pm,B) “t ^Pm,BPm+l,Bi,^ Pm,A) — O 5 

QPm,Ai^ PmAl,A) QPm+l,A{^ Pm+2,A) “t“ UJdi^K {^K -j- 

^Pm+l,APm+ 2 ,A (1 Pm+1,B) “t ^ Pm+1,B Pm+2,B{)^ Pm+1,A^ — O 5 

2 J^2 ~ — {K + 1)pb) + ^^Pa(1 “ Pb) — ^P%{^ — Pa) = 0 , 

where 0 < S' < me, {m + l)e < x < 1 and 0 < x < 1. 

Since it is difficult to solve the system (|9]) analytically due to its hybrid 
nature, we propose to hnd the steady-state solution of the hybrid system 
using the singular perturbation technique Recently, this technique has 
been successfully applied for studying two-channel homogeneous TASEP with 
LK 0, 10. 

Singular perturbation technique typically involves obtaining solutions of 
the differential equations describing the boundary layer region and the bulk 
solution, separately and then matching both the solutions to get the global 
solution. The bulk part of the solution, known as outer solution, is found 
in the limit e —>■ 0. We propose a generalized methodology to hnd the 
outer solution of the hybrid system ([9]). Instead of solving the system ([9]) 
explicitly, we hnd the steady-state outer solution by capturing th& Ipsg time 
solution of ([6]), ([7]) and (jH]) using the scheme discussed below [6, 


lOl. We 


discretize the continuum part of model equations using hnite-diherence and 
replace the time derivative in discrete equations at sites (m. A) and A) 

with forward-diherence formula. For l<i<m, m-l-l<i<L with 
j = A and 1 < i < L with j = B, the following scheme is applied to model 
equations o, o and dH]) with the boundary conditions pa(0) = Pb(0) = a 
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and pa( 1) = Ps(l) = 1-13. 





1 


( 10 ) 


Pm 




The solution is captured in the limit n —>■ cxo to ensure the occurrence of 
a steady-state. Clearly, the above system is consistent and gives a unique 
steady-state outer solution. It is easy to check the consistency of the above 
system of equations by analyzing the system in the limit g —)■ 1. 

To satisfy the boundary conditions at both the ends, the density prohles 
incurs a crossover narrow regime in the form of either a boundary layer or 
a shock. This solution is known as inner solution and is found by rescaling 
the space variable in the neighbourhood of Xd as x = where Xd is the 

position of boundary layer. This rescaling eliminates the non-conservative 
source terms, which are formed by lane-changing transitions and attachment- 
detachment dynamics, in the system ([2]). In terms of x, the inner solution 
Pj^in is given by 



( 11 ) 


Here, the integration constant aj is computed from the matching condition 
of outer and inner solution. For example, let the boundary layer appears at 
right boundary (a; = 1) in lane j, the matching condition requires pj^inix —>■ 
—oo) = pj^out{x = 1) = Pj^o- Here, pj^o is value of the outer solution in lane-j 
at a; = 1. Clearly, pj^o is a function of system parameters Qd, and 
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Solving Eq. ([TT]) with aj = — pj^o] we get 


ft,in = i — !i tanh I x|2p^,„ - 1| + 0 ;. 


( 12 ) 


where = tanh 


-1 


1 - 2,5 


The solution given by eq. flT^ represents a 


y|2pj,o 1| 

right boundary layer (rbl) in lane j with positive slope (tanh—r). When 
/3 < pj^oi the inner solution fails to satisfy the right boundary condition 
Pj,in{x —)■ cx)) = 1-/3 and deconhnes from the boundary to enter the bulk 
of lane-j in the form of a shock. Thus (3 = Pj^o{oi) acts as a bulk phase 
transition line. Apart from deconhnement, the inner solution also undergoes 
change in its slope across the line /3 = 1 — pj ^(a). The right boundary layer 
with negative slope is given by pj^in = | + coth (^\2pj^o “ 1 | + , 




1 - 2/3 


Here, the change in the slope of boundary 


where = coth , ,75 -^ 

layer describes a surface transition, which does not affect bulk density prohle. 
Analyzing, in a similar fashion, we can End mathematically all the bulk as 
well as surface transition lines in the a — f3 plane and generate the phase 
diagram for different parameters. For a more detailed solution methodology 
of singular perturbation technique, we refer 6 |. 


5. Monte Carlo simulations 

To testify the outcomes of our theoretical mean-field approach, we sim¬ 
ulate the two-channel model using Monte Carlo simulations (MCS) jlj]. We 
adopt random sequential update dynamical rules, which is the realization of 
the usual master equation in continuous time. Moreover, random-sequential 
type update rule is appropriate to model the dynamics of intracellular protein 


transport 


25 


, l 6 |. We run the Monte Carlo simulations for 10^° — 10^^ time 
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steps and first 5% steps are ignored to ensure the occurrence of steady-state. 
The densities in both the lanes have been computed by taking time averages 
over an interval of lOL. Since the size of a real system is normally not very 
large, it is reasonable to simulate the system for a lattice size up to L = 1000. 
Note that the dynamical rules involve many kinds of particle hoppings such 
as horizontal hoppings, vertical transitions and attachment-detachment dy¬ 
namics, it is computationally expensive to calculate steady-state densities. 
As we proceed in the next section, it is found that the results of theoretical 
mean-held approach agree to a good extent with Monte Carlo simulations. 


6. Results and discussion 


In this section, we provide a detailed discussion of the steady-state prop¬ 
erties of the proposed model, with emphasis on the resulting phase diagrams 
and analyze the effects of bottleneck strength, lane-changing rate and binding 
constant. For simplicity, we assume that the bottleneck in lane A is located 
exactly at the middle i.e. m = L/2. The steady-state phase diagrams as 
well as density prohles are obtained from hybrid mean-held approach and 
validated using extensive Monte Carlo simulations. 

When there is no bottleneck, i.e. q = 1, the present model reduces to 
the homogeneous two-channel TASEP with LK. Under symmetric coupling 
conditions, the two-channel homogeneous TASEP with LK is analogous to its 
single-channel counterpart, which has been studied well in literature 3^. For 
g > 1, we call the defect to be fast defect as the bottleneck rate is more than 
the usual horizontal particle hopping rate in the system. Clearly, a single fast 
defect does not ahect the steady-state dynamics as the rate-limiting hoppings 
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occurring with rate 1 dominate the whole lattice 18| . Thus, it is appropriate 
to focus on the situation where the bottleneck acts as a slow defect i.e. q <1. 

It is well known that the steady-state dynamics of exclusion process are 
signihcantly affected by the binding constant. So, we proceed by segmenting 
the analysis into two subcases viz. (i) K = 1, and (ii) K ^ 1. The other 
parameters are kept hxed to VLd = 0.2 and L = 1000. 

6.1. K = 1 

We examine the steady-state dynamics for different strengths of bottle¬ 
neck by gradually decreasing the value of q. No signihcant qualitative changes 
are found in the phase diagram till gc,i ~ 0.75, below which one notices the 
effect of bottleneck more clearly. Fig. |2] shows the phase diagram of the 
system with q = 0.75 and = 1. Specihcally, £g. Ha) shows the complete 
characterization of various phases in terms of the boundary layers. The solid 
and dashed lines denote the bulk and surface transitions, respectively. As an 
illustration, the notation (tanh-r,tanh-r) signifies that there is a right bound¬ 
ary layer in lane A as well as lane B of the kind tanh-r, given by eq. flT^ and 
S denotes the shock phase. For a detailed understanding of the bulk and 


surface transitions, we refer the readers to references 


23| . As discussed 


in previous section that a surface transition is responsible for the change in 
slope of the boundary layer only and does not affect the bulk density, hence¬ 
forth, our focus will be only on bulk transitions to avoid any complexities 
in the phase diagrams. Fig. represents the bulk phase transitions only 
corresponding to £g. Ha). 

To understand the dynamics in a better way, we introduce an appropriate 
notation {PhA{ieft) — PhA(right), P^b) denoting the kind of stationary phase 


16 






(a) 




Figure 2: For q = 0.75, ^Id = 0.2, 12 = 1 and L = 1000, (a) Phase 

diagram showing bulk as well as surface transitions and classification accord¬ 
ing to the boundary layers, NBL denotes no boundary layer, (b) Phase dia¬ 
gram showing qualitative nature of phases with no surface transitions. Notations:- 
M: Meissner phase, 1: (S'6-LD/MC/HD,LD/HD/LD/MC/HD), 2: (LD/MC/HD- 
LD/MC/HD,LD/MC/HD/LD/MC/HD), 3: (LD/MC/HD-S,LD/MC/HD/LD/HD), 

4: (S'b-S,LD/HD/LD/HD), 5: (S'b-LD/MC,LD/HD/LD/MC), 6: (LD/MC/HD- 

LD/MC,LD/MC/HD/LD/MC), 7:(MC/HD-LD/MC,MC/HD/LD/MC), 8: (MC/HD- 
LD/MC/HD,MC/HD/LD/MC/HD) and 9: (MC/HD-S,MC/HD/LD/HD). Solid and 
dashed lines denote bulk and surface transitions, respectively. 


in the two-channel system. Here, PhA(ieft){.PhA{right)) denotes the type of 
stationary phase in left (right) subsystem of lane A; while Phs denotes the 
stationary phase in lane B. For example, the phase (LD-S,S) denotes that 
the region left to the bottleneck is low density (LD) with a shock present in 
the right subsystem in lane A and a shock is also present in lane B. 

We list important features of the steady-state dynamics under the current 
choice of parameters. 

1. Meissner phase and bottleneck affected region: In the upper right quad- 
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(b) 





Figure 3: Density profiles for q = 0.75, fid = 0.2, D = 1 and L = 1000. (a) (LD- 
LD,LD/LD) phase for a = 0.2,/? = 0.8, (b) (LD-S,LD/S) phase for a = 0.1,/3 = 0.2, 
(c)(HD-HD,HD/HD) for a = 0.5,/3 = 0.15, (d) (S-HD,S/HD) phase for a = 0.2,^ = 
0.1, (e) (MC/HD-LD/MC,MC/HD/LD/MC) phase for a = 0.8,/3 = 0.8 and (f) 
S,LD/HD/LD/HD) phase for a = 0.4,^ = 0.4. 
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rant (a > 0.5,/? > 0.5, shaded blue), the bulk densities in each lane 
becomes independent of entrance as well as exit rates. This effect is 
similar to the behavior of the Meissner phase found in superconducting 
materials 12|]. Therefore, we also adopt the name Meissner phase (de¬ 


noted by M) to identify boundary independent behavior of two-channel 
system. It is important to mention that Meissner phase is actually a 
part of certain existing phase, denoted by phase 7 in £g. |2](b). Meissner 
phase is also reported in the literature for a homogeneous single-channel 


TASEP with LK 


25|. Note that the density prohles get signihcantly 


influenced by bottleneck for a specihc range of a, /5. The corresponding 
region in phase-plane is marked with red dots and referred as bottle¬ 
neck phase throughout this paper. Similar results also exist for single¬ 
channel inhomogeneous TASEP with LK, in which bottleneck phase is 
identihed and analyzed in terms of carrying capacity 2^. In contrast, 
for low entrance and exit rates one recovers the TASEP/LK density 
prohles perturbed by a local spike or a dip, which will be discussed as 
we proceed in this paper. 

2. Phase transitions and mixed phases: To gain an intuitive understand¬ 
ing of the steady-state dynamics of the system, it is important to an¬ 
alyze the phase transitions. Across the bulk transition line between 
(LD-LD,LD)(fig. |3](a)) and (LD-S,S) (kg. [3](b)) phases, one finds the 
formation of a shock through deconfinement of the right boundary 


layer 


10|. It is evident from £g. |2](b) that (d = a + ^ld = o: + 0.2 


and (3 = a — VLd = (y. — 0.2 for a, /9 < 1/2 are the lines of deconhnement 
of shock from right boundary layer (LD phase) and left boundary layer 
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(HD phase), respectively. Interestingly, these are also the correspond¬ 
ing lines of deconhnement in analogous single-channel homogeneous 


system 2^. From here, one can infer that the phase diagram under a 
weak bottleneck has almost similar topological structure as for the one 
with no bottleneck {q = 1) 


25 


though some of the phases differ 


in their qualitative nature. Actually, a weak bottleneck is incapable 
of affecting the steady-state dynamics of a symmetrically coupled two- 
channel TASEP with LK system globally. The only effects can be seen 
in the vicinity of the bottleneck. This point is evident from the density 
prohles in Eg. [3l which will be discussed in detail subsequently. 

As one can see from Eg. [3](b), a shock is present in both the lanes in 
(LD-S,S) phase. On decreasing (3, the shock in both the lanes moves 
leftwards and reaches exactly at x = 1/2 on the line a = /3; a, 13 < 1/2. 
This line corresponds to the phase transition between (LD-S,S) and (S- 
HD,S) phases. A further reduction in magnitude of (3 moves the shock 
upstream to the bottleneck giving rise to (S-HD,S) phase. The bulk 
phase transition from (HD-HD,HD)(Eg. |3](c))to (S-HD,S)(Eg. |3](d)) can 
also be understood on the similar lines. Rest of the phase transitions 
in Eg. [2](b) are due to gradual changes in the magnitude of densities 
with respect to entrance and exit rates. 

Note that lane B has exactly the same density as that in lane A every¬ 
where, except in the vicinity of the bottleneck (see Eg. |3]). Interestingly, 
this slight difierence in density near the bottleneck is able to produce 
peculiar kinds of density proEles, giving rise to a number of coexistence 
(or mixed) phases such as (LD/MC/HD-S,LD/MC/HD/LD/HD) and 
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(MC/HD-LD/MC,MC/HD/LD/MC) etc. An illustration of a mixed 
phase is shown in hg. |3]^e). 

3. Formation of upward and downward spikes: Consider the LD phase 
in left subsystem of lane A, where the particles have enough space to 
move along the lattice and the presence of a weak bottleneck does not 
affect the particle motion significantly. Here, the density profile in lane 
A incurs a local perturbation in the form of an upward spike as shown 
in fig. El^a). Due to comparatively slower particle hopping rate at the 
bottleneck, the particles coming from the left feels a little congestion 
before entering the bottleneck. Due to this, the density in lane A 
rises slightly just upstream to the bottleneck, leading to formation of 
an upward spike in the density prohle. Importantly, lane B remains 
unaffected in this phase. Similarly, one can understand the formation 
of a downward spike shown in fig. [3]^b)) in HD phase. This can be 
attributed to particle-hole symmetry. When a particle comes out of 
the bottleneck in a highly packed situation, it experiences a sudden 
increase in its hopping rate. This leads to formation of a lower density 


region just downstream to t 
reported in literature j3. 


37| 


re bottleneck. Similar hndings have been 


4. Transition to bottleneck-induced shock: Now we discuss an important 
feature of our system dynamics i.e. the formation of a bottleneck- 
induced shock. With LD in lane A, an increase in entrance rate a 
increases the number of particles in left subsystem of lane A, which 
ultimately increases height of the existing upward spike. At a par¬ 
ticular a = ttc = 0.3, this local density perturbation converts into a 
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shock, which travels in the left subsystem in lane A for further increase 
in a. Since this shock emerges due to bottleneck only, we call this 
shock as bottleneck-induced shock and denote it by Sb- The formation 
of bottleneck-induced shock is also observed for a single-channel inho¬ 
mogeneous TASEP with LK . The presence of bottleneck-induced 
shock leads to the existence of three new phases (S'b-S,LD/HD/LD/HD), 
(S'fe-LD/MC/HD,LD/HD/LD/MC) and (Afe-LD/MC/HD,LD/HD/LD/ 
MC/HD) in the phase-plane as shown in Eg. [2]^b). The existence of 
mixed phases in lane B indicates the effect of bottleneck in homoge¬ 
neous lane as well. Interestingly, there exist two shocks in lane A in {Sb- 
S,LD/HD/LD/HD) phase (See £g. |3](f)), the left one being bottleneck- 
induced; whereas the right one is formed due to deconhnement of right 
boundary layer. 

6.1.1. Effect of bottleneck strength 

We propose a terminology to identify the strength of the bottleneck as 
shown in table [H We also provide the number of steady-state phases observed 
under each case, the details of which follow in the text. 

Weak bottleneck: So far, we have discussed the steady-state dynamics for 
a weak bottleneck and pointed out important features such as appearance of 
Meissner phase, bottleneck phase and bottleneck-induced shock etc. Clearly, 
the number of steady-state phases under a weak bottleneck is more than in 
homogenous case, which is due to the appearance of mixed phases. Now, we 
will gradually change the strength of bottleneck from weak to moderate and 
then to strong and see whether the same trend continues or not? 

Moderate strength: In the range of moderate strength of bottleneck, we 
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Table 1: Steady-state properties of the system with fid = 0.2, = 1. if = 3 is chosen as 

a specific value for the case K ^ 1. 


Strength 

Bottleneck 

Number of stationary phases 


rate 

K 1 

K = 1 

No bottleneck 

q = l 

4 

7 

Weak 

0.75 < g < 1 

2 

13 

Moderate 

0.4 < g < 0.75 

2 

9 


0.1 < g < 0.4 

4 


Strong 



14 


0 < g < 0.1 

3 


Blockage 

g = 0 

3 

14 


have chosen q = 0.5 to generate the phase diagram. It is found that for 
a moderate bottleneck, the phase diagram (See £g. 111(a)) is qualitatively as 
well as quantitatively different from the one with weaker bottleneck (See 
£g. [2](b)). An important distinguishing feature of the phase diagram is the 
appearance of bottleneck-induced shock in lane B as well. This can be seen 
from £g. 111(a), which identifies (S';,-S,S'b/S) as a new phase. While the nature 
of bottleneck-induced shock in lane A has already been discussed, one needs 
to investigate the occurrence of two shocks in lane B. Clearly one located in 
the right is due to the deconhnement of right boundary layer; while other is 
because of the indirect impact of bottleneck in lane A, whose origin can be 
understood as follows. Actually, the presence of a bottleneck-induced shock 
in lane A produces a local high-density region there, making it difficult to 
accommodate incoming particles from lane B. As a result, the density in the 
left subsystem of lane B rises and converts into a shock with an increase in 


23 










entrance rate a. The shock produced in lane B, due to the effect of bottleneck 
in lane A, is also called as bottleneck-induced shock. 



Figure 4: Phase diagrams for fid = 0.2, ft = 1 and L — 1000 with {a.)q = 0.5 
and (b)density profile in {Sb-S,Sb/S) phase with a = 0.4, /3 = 0.4 and q = 0.5, and (c) 
phase diagram for q = 0.1, I: {Sb-S,Sb/S), II: (S'b-LD,S'(,/S), III: (HD-LD,S'h/S) and IV: 
(HD-S,5'b/S). Increase in height of upward spike (shown in magnified inset) and conversion 
of spike into bottleneck-induced shock with a decrease in q 


On moving from weak to moderate bottleneck, the phase diagram gets 
simplihed resulting into a reduction in number of phases. The coexistence 
phases (MC/HD-LD/MC,MC/HD/LD/MC), (MC/HD-S,MC/HD/LD/HD) 
and (S'fc-LD/MC,LD/HD/LD/MC), observed in £g. [2](b), are replaced by 
(HD-LD,HD/LD), (HD-S,S) and (S'fe-LD,S), respectively. This shows that 
MC phase, either as a whole or as a part of density profile, does not exist 
for moderate bottleneck, contrary to the weak bottleneck case. This result is 


analogous to t 


bottleneck 
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le one found for single-channel TASEP with LK with a single 
Further, the bottleneck phase (dotted region) enlarges, as 


expected for a stronger bottleneck; whereas the phases (LD-LD,LD),(LD- 
S,S), (S-HD,S) and (HD-HD,HD) shrink. Meissner phase remains unaffected 
with increase in bottleneck strength. 
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Strong bottleneck: Fig. SJ^c) shows the phase diagram for the present 
system with q = 0.1, a strong bottleneck. It can be seen that topology of 
the phase diagram is signihcantly changed in comparison to the one with 
moderate bottleneck. Here, the bottleneck effect is sufficiently strong to 
increase the density difference between two lanes up to a level that we have 
a complex phase diagram structure as shown in £g. ID^c). One can see that 
LD phase does not exist in the left subsystem to the bottleneck in lane A. 
This is physically justihed as the extremely slow hopping rate at bottleneck 
obstructs the incoming particles, leading to formation of either a shock or a 
HD region in left subsystem. Similar argument and particle-hole symmetry 
accounts for the non-existence of HD phase downstream to the bottleneck in 
lane A. It further accounts for the enlargement of (HD-LD,HD/LD) phase 
with a decrease in q. Note that the bottleneck phase (dotted region) further 
enlarges to maximum extent and covers the whole a — /3 plane. 

We have also checked the results for the case of complete blockage at 
the bottleneck i.e. q = 0. The phase diagram for q = 0 remains essentially 
similar to the one for q = 0.1. The only change is shifting of the phase 
boundaries, which leads to expansion of (HD-LD, HD/LD) phase (following 
the similar trend as explained above) and shrinkage of rest of the phases. 

Fig. |5] shows the transition of upward spike into bottleneck-induced shock 
with respect to an increase in bottleneck strength. Firstly the height of up¬ 
ward spike in lane A at a: = 1/2 increases in moderate bottleneck case(clearly 
seen in the inset), then the spike converts into shock for strong bottleneck. 
A further decrease in value of q pulls the shock leftwards. 
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Figure 5: Increase in height of upward spike (shown in magnified inset) and conversion of 
spike into bottleneck-induced shock with a decrease in q. 


6.1.2. Effect of coupling strength 

So far, we have investigated the system dynamics under different bottle¬ 
neck strengths, but for hxed lane-changing rate at = 1. Since coupling 
strength has always been an important parameter in two-channel systems, it 
is necessary to examine its role for the proposed model. In literature, it is 


well explored 


al9 


28 


36| that higher orders of lane-changing rate does not 


affect the dynamics of a symmetrically coupled homogeneous two-channel 
TASEP with/without LK because the system decouples into two indepen¬ 
dent single-channel homogeneous systems. But, in our model, this is not 
the case as our system is inhomogeneous. So, we thoroughly investigate the 
effect of coupling strength on the dynamics of the proposed inhomogeneous 
model under different strengths of bottleneck systematically. Our aim is to 
discuss the important topological changes in the phase diagram induced by 
higher order of lane-changing rate (For dehnition of order, please see ^). 

For a weak bottleneck, the effect of lane-changing rate on the phase di- 
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agram is insignificant till (9(f2) < 1000(^2^). Only little qnantitative differ¬ 
ences are seen on density profiles in both the lanes. When 0(0) > 1000(0^), 
we notice signihcant changes in the density profiles of both the lanes leading 
to some qualitative changes in the nature of certain phases. Fig. [6](a) shows 
the phase diagram for q = 0.75 with O = 100. For higher orders of O, we have 
found that the number of stationary phases remain same. Interestingly, some 
of the existing phases get transformed into new phases without a noticeable 
change in the phase boundaries. For example, the phases (LD/MC/HD- 
LD/MC,LD/MC/HD/LD/MC), (MC/HD-LD/MC,MC/HD/LD/MC) and 

(5b-S,LD/HD/LD/HD) change to (LD/HD-MC,LD/HD/MC), (MC/HD- 
MC, MC/HD/MC), (S'b-HD,5'6/HD), respectively on moving from hi = 1 to 

n = 100 . 

The effect of increasing the order of 12 on a moderate bottleneck is sim¬ 
ilar to the one discussed for weak bottleneck. Therefore, we skip the phase 
diagrams for this case at higher orders of 12. 

Moving our attention towards strong bottleneck, we analyze the phase 
diagram for g = 0.1. No structural changes in the phase diagram are seen 
till 0(12) < 10O(12rf), except minor translation of various phase boundaries. 
As an example, the (S'fe-S,S) phase expands; while (HD-LD,HD/LD) phase 
shrinks. An important phenomenon occurs at 0(12) = 100(12^), where the 
density difference between both the lanes, in the vicinity of the bottleneck, de¬ 
creases. As a result, the complexity of the phase diagram reduces as shown in 
fig.[6](b). Importantly, four new phases (LD-LD,LD), (LD-S,S), (S-IID,S) and 
(HD-HD,HD) appear. The effect of bottleneck, despite its higher strength at 
q = 0.1, reduces because the bottleneck phase covers comparatively smaller 
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region in a—f3 plane. Further increase in lane-changing rate does not produce 
any subsequent changes in the phase diagram for q = 0.1. This happens due 
to the weakening of bottleneck effect by an increase in lane-changing rate, 
which is discussed as follows. 
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Figure 6: Phase diagram for 12^ = 0.2, and L = 1000 with (a) q = 

0.75, n = 100, Notations:- 1: (S'6-LD/MC/HD,S'6/LD/MC/HD), 2: (LD/HD- 

LD/MC/HD,LD/HD/LD/MC/HD), 3: (LD/HD-HD,LD/HD) and 4: (S'6-HD,S'b/HD).(b) 
g = 0.1 with f2 = 10. 


Weakening of bottleneck effect It can be seen from hgs. [21 0] and [6] that 
the bottleneck affected region (dotted region) shrinks at higher orders of 
lane-changing rate in comparison to the one for lower orders of under all 
values of q. Specihcally, the weakening of bottleneck effect occurs for lower 
values of a and ft with an increase in lane-changing rate fl. This is further 
justihed from £g. [Tj, which shows a monotonic reduction in peak’s height of 
upward spike at bottleneck with respect to an increase in ff. An obvious 
explanation of this effect is that the particles in lane A, feeling hindered due 
to the bottleneck, would shift to other lane more frequently. As a result, the 
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congestion in lane A gets slightly reduced as one increase the 
rate. Similar kind of effect is also reported in literature 


lane-changing 
The insets in 


hg. [7] shows the similar effect for downward spike as well. 



0 200 400 600 800 1000 
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Figure 7: Variation in height of upward spike (insets show for downward spike) at bottle¬ 
neck in lane A for (a) q = 0.75, (b) q = 0.5 and (c) q = 0.1. Here, a = 0.1, /3 = 0.6 for 
upward spike in LD phase and a = 0.8, /3 = 0.1 for downward spike in HD phase. 


6.S. K ^ 1 

Since the methodology of the proposed hybrid mean-held approach is 
general and applicable for any value of binding constant K, we follow same 
lines to discuss the case where attachment and detachment rates differ, ^ 
ild, i-e. K 1. 

Fig. [8] shows the steady-state phase diagrams with K = 3. Comparing 
hg. [8]with hgs. [2]and HI we notice that the complexity of the phase diagrams 
reduces tremendously for iC = 3 under all values of q. The number of steady- 
state phases is shown in table [H Note that due to the increased attachment 
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(a) 
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Figure 8: The phase diagrams for K = 3, fid = 0.2 and = 0.5 with different q: (a) 
q = 0.75, (b) q = 0.4, and (c) q = 0.1. Notations, M: Meissner phase (shaded blue) and 
BP; Bottleneck phase (shaded with red dots) 


rate in comparison to detachment rate, the density in both the lanes increases, 
leading to a complete absence of LD phase in the system. 

For weaker bottleneck, the phase diagram is much simple consisting of 
only two stationary phases (S-HD,S) and (HD-HD,HD) as shown in £g. [Hl^a). 
Quite surprisingly, no density prohle in a — /? plane is signihcantly affected by 
bottleneck, excluding the presence of a downward spike at a; = 1/2. This is in 
sharp contrast to the parallel case with K = 1 (See £g. [2])b)). Consequently, 
one does not get a bottleneck-induced shock here for any values of a and /3. 

Further decrease in the value of q does not change the structure of phase 
diagram but increases the height of downward spike in HD phase. This trend 
continues till one reaches gc ,2 ~ 0.4, at which the perturbation of downward 
spike enters the right subsystem in lane A in the form of a shock, called as 
bottleneck-induced shock. Its emergence leads to formation of new steady- 
state phases such as {S — Sb, S) and (HD-S'f„HD) as shown in £g. IH](b). A 
part of phase-plane (shaded with red dots) is strongly influenced by the 
bottleneck-induced shock called as bottleneck phase, which was absent for 
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q > 0.4 (See fig. [Hl(a))- On the other hand, the value of q does not influence 
the Meissner phase. 

A further decrease in value of q below gc ,2 leads to signihcant changes 
in the steady-state dynamics. The emergence of a new phase (HD-S'6,S) is 
seen at another critical value gc ,3 = 0.1 as shown in £g. [Hl(c). Additionally 
two of the earlier existing phases (S-HD,S) and (HD-HD,HD) disappear com¬ 
pletely; while the bottleneck phase covers the whole a — /3 plane reflecting 
the profound effect of strong bottleneck. 

6.2.1. Effect of coupling strength 

We have also investigated the effect of lane-changing rates on the dynam¬ 
ics of the proposed system under case K 1. For 0(fl) < 100(11^), no 
signihcant changes are seen except shifting of phase boundaries for any value 
of q. We focus on the situation when 0(f2) > 1000(11^). It is found that 
the steady-state phase diagram for weak bottlenecks remains unchanged and 
has a structure similar to the one seen in £g. [Hl(a)- But the height of down¬ 
ward spike at x = 1/2 decreases with respect to increase in lane-changing 
rate, analogously to the parallel case K = 1. One can infer from here that 
an increase in lane-changing rate weakens the bottleneck effect with K ^ 1 
as well. The extent of weakening effect is up to a level that we do not get 
a bottleneck-induced shock till q > 0.15, which means the global effects of 
the disturbance caused by bottleneck can be captured only for a very strong 
bottleneck. As an illustration, we provide the phase diagram for g = 0.1 in 
£g. [n](a), which consists of four new stationary phases. Moreover, the bottle¬ 
neck phase shrinks in its size, which once again indicates the weakening of 
the bottleneck’s effect with respect to an increases in fl. 
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Figure 9: Phase diagrams for (a) q = 0.1, K = 3, = 100 and {h)q = 0.01, K = A, 

n = 1000 . 


6.3. Non-existence of bottleneck-induced shock 

We have observed that an increase in lane-changing rate weakens the 
bottleneck effect and as a result the bottleneck affected region m a — (3 
plane shrinks with respect to an increase in kl. Comparing the two cases 
K = 1 and K ^ 1, one can infer that the effect of bottleneck is further 
reduced when K ^ 1. The above two statements generate interest to know 
whether there exist some values of K and/or at which we do not get a 
bottleneck-induced shock even for a strong bottleneck. A thorough analysis 
of our study has revealed that for K = 4 and k2 = 1000 with a very strong 
bottleneck (g = 0.01, almost a blockage), we get a steady-state situation free 
of any bottleneck-induced shock. The resulting steady-state phase diagram 
is shown in fig. |9])b). Although we do have a downward spike in the density 
profiles, yet the non-existence of a bottleneck-induced shock is an important 
outcome. Such dynamics further opens up new challenges to dig deeply into 
the shock dynamics of the proposed model, which will be undertaken in our 
future study. 
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7. Conclusion 


In this paper, we attempted to provide a clear picture of the role played 
by an inhomogeneity in a two-channel TASEP with LK under a symmet¬ 
ric coupling environment. The location of bottleneck is hxed at middle site 
in lane A, while lane B is kept homogeneous. Theoretically, we have in¬ 
troduced a new hybrid approach and also employed singular perturbation 
technique to get the steady-state density prohles and phase diagrams. We 
have adopted random-sequential update rules to perform Monte Carlo simu¬ 
lations and shown a good agreement between theoretical and MCS. 

The study is segmented into two subcases: K = 1 and K 1. Based 
on the analysis, one finds very interesting topologies of the phase diagrams 
under different bottleneck rates. It is concluded that the steady-state phase 
diagrams for different bottleneck strengths under K = 1 are comparatively 
richer in their composition as compared to those under K ^ 1. For lower 
strengths of bottleneck, the phase diagram comes out to be quite similar 
to the corresponding homogeneous system and effects on the density pro¬ 
hles are only local in the vicinity of the bottleneck. The origin of upward 
and downward spike at bottleneck in LD and HD phases, respectively is ex¬ 
plained. We also shed light on the emergence of bottleneck-induced shock in 
the system, which accounts for the global effects of bottleneck. The effect of 
lane-changing rate on the dynamics is also examined and found that an in¬ 
crease in coupling strength weakens the bottleneck effect for any value of K. 
This effect is supported by the reduction in height of upward and downward 
spikes with respect to increase in lane-changing rate. Another important ob¬ 
servation is that unequal attachment-detachment rates are also responsible 
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for weakening of bottleneck effect. 

Here, we focused on a localized bottleneck, fixed at the middle of the 
lattice far away from the boundaries. The approach developed in this paper 
might be easily extended to analyze more general inhomogeneous systems 
such as those with more number of bottlenecks and also to study the influ¬ 
ence of the location of the bottleneck in the multi-channel TASEPs. From 
the biological signihcance, one can say that the present study might help 
in optimizing the protein translation rate in multi-channel inhomogeneous 
systems. 
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